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Abstract 

We perform a complete one-loop computation of the two-body flavour-changing top 
decays f —)■ ch and f —)■ cU (V = j, Z), within the aligned two-Higgs-doublet model. 
We evaluate the impact of the model parameters on the associated branching ratios, 
taking into account constraints from flavour data and measurements of the Higgs 
properties. Assuming that the 125 GeV Higgs corresponds to the lightest CP-even 
scalar of the CP-conserving aligned two-Higgs-doublet model, we hud that the rates 
for such flavour-changing top decays he below the expected sensitivity of the future 
high-luminosity phase of the LHC. Measurements of the Higgs signal strength in 
the di-photon channel are found to play an important role in limiting the size of the 
t —)■ ch decay rate when the charged scalar of the model is light. 



1 Introduction 


The discovery of a new boson, with a mass close to 125 GeV, by the ATLAS and CMS 
collaborations stands as a remarkable success of the Standard Model (SM) of electroweak in¬ 
teractions. The properties of this boson are so far in agreement with those of the SM Higgs, 
indicating that the new particle is indeed associated with the mechanism of electroweak sym¬ 
metry breaking (EWSB). A characteristic feature of the SM is the absence of flavour-changing 
neutral-current (FCNC) interactions at tree-level. FCNCs are generated through quantum loop 
corrections in the SM, but they are strongly suppressed by the Glashow-Iliopoulos-Maiani 
(GIM) mechanism [^. 

A widely studied enlargement of the electroweak theory consists in adding a second scalar 
doublet to the SM held content. The so-called two-Higgs-doublet model (2HDM) represents 
a minimal extension of the SM scalar sector that easily accommodates electroweak precision 
data and leads to a very rich phenomenology [^. In the most general version of the 2HDM 
unwanted FGNGs appear at tree-level, which represents a major shortcoming of the model. The 
hypothesis of natural havour conservation (NFG) is the usual way out to this issue. By limiting 
the number of scalar doublets coupling to a given type of right-handed fermion to be at most 
one, the absence of dangerous FGNGs is guaranteed [^[^. A more general solution is that of 
Yukawa alignment (^. The aligned two-Higgs-doublet model (A2HDM) assumes that the two 
Yukawa matrices coupled to the same type of right-handed fermion are aligned in flavour space, 
so that no FGNGs appear at tree level. Explicit models where a Yukawa aligned structure arises 
due to an underlying symmetry have been discussed in refs. [8 -12 . Interestingly, all different 


versions of the 2HDM with NFG are recovered as particular limits of the A2HDM. Gonstraints 
on the A2HDM from flavour and collider data have been analyzed in refs. 13-18 and 19 26 , 
respectively, extracting relevant bounds on the model parameters. 

In this work we study the flavour-changing top-quark decays t ^ ch and t ^ cV (V = 7 , Z), 
within the framework of the A2HDM. They arise at the loop level and are strongly suppressed in 
the SM, due to the GIM mechanism. A signihcant enhancement can be achieved in alternative 
scenarios of EWSB, making these processes a suitable place to look for new physics beyond 
the SM. Early considerations of these effects were done in refs. 27-32 . A concise review of 
the flavour-changing top-decay phenomenology can be found in ref. 33 . In the A2HDM these 


decays receive additional charged Higgs contributions at the one-loop level which could lift the 
decay rates. 

Gomprehensive analyses of flavour-changing top decays within 2HDMs with NFG have been 
done in refs. |^[3l|[Ml|^ . However, these rare processes have not been investigated yet within 
the more general setting of the A2HDM. Furthermore, since these studies were performed, 
considerable experimental progress in our understanding of the EWSB mechanism has been 
made, translating into tight constraints on possible extensions of the SM scalar sector. 


Searches for flavour-changing top decays have been performed recently by the ATLAS 36 37 


and GMS collaborations 38 -40 , placing limits on the associated branching ratios. An overview 
of the current experimental status can be found in ref. 41 . With the large amount of data 


that will be collected in the future LHG runs, it is expected that these bounds will be improved 
by at least one order of magnitude [42,43 . 


Our paper is organized as follows. In section we introduce the A2HDM. Flavour-changing 
top decays are discussed in section [^ which presents the results of our calculations. A phe- 
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nomenological analysis of these processes is given in section and onr conclnsions are finally 
snmmarized in section Explicit analytical resnlts for the relevant decay amplitndes are given 
in the appendices. 


2 Framework 


The 2HDM extends the SM scalar sector with an additional complex scalar donblet. In the 
Higgs basis, where only one donblet acqnires vacnnm expectation valne, the scalar helds are 
parametrized by 


21 


$1 = 


G'+ 

-j,{v + S, + iG°) 


L V2 


$2 — 


H+ 

(S 2 + 2 S:,} 


( 2 . 1 ) 


with V = ~ 246 GeV. Here correspond to the wonld-be Goldstone bosons, 

giving mass to the gange vector bosons, while is a charged Higgs. The scalar spectrnm 
also contains three nentral Higgs bosons = {h{x), H(x), A{x)}, given by (p° = TZjkSk, 

where TZ is an orthogonal matrix obtained after diagonalizing the mass terms in the scalar 
potential [^. In general none of the nentral Higgs bosons are GP eigenstates. 


2.1 Scalar sector 


The most general scalar potential allowed by the electroweak gange symmetry can be written 
as 


V = III + p 2 ^1^2 + h3 ^[^2 + hs $1 


Ai + As <hT<f )2 + As ^^2 + A 4 <hl4>2 


A 5 T Ag + A 7 $ 2 *^ 2 ^ ^^ 1 *^ 2 ^ + h.c. 


( 2 . 2 ) 


Dne to the Hermiticity of the scalar potential, all parameters are real with the exception of 
hs; As, Ag and A 7 . The minimization conditions impose the relations 


/ii — —Ai , p-3 — 

The mass of the charged Higgs is given by 


— ps + 2 A3 ) 



(2.3) 


(2.4) 


while those of the nentral scalars have been obtained in ref. 21 to hrst order in the GP-violating 
parameters. In the GP-conserving limit the nentral Higgs bosons are GP eigenstates. The GP- 
odd held A corresponds to S 3 and the GP-even states are orthogonal combinations of Si and 


S 2 -. 



COS a sin a 
— sin a cos a 



(2.5) 
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( 2 . 6 ) 


Here ^ hy convention and the mixing angle a is determined by 


sin 2 a = 


— 2 Xqv‘^ 

Mh-Mr 


cos 2(5 = 


Ml+ 2 (A 5 -Aiy 
Mh-Ml 


In the CP-conserving limit the masses of the neutral Higgs bosons are given by 

M,; = i(E-A), M==i(E + A), AfJ = Aj ) , 


where 


E = MU + I 2Ai + h + As I , 


A = 




(2.7) 


( 2 . 8 ) 


(2.9) 


By performing a phase redehnition of the CP-even helds one can restrict the mixing angle to 
the range 0 ^ d < vr. Moreover, the Higgs basis of the CP-conserving 2HDM is dehned up to 


a global rephasing of the second Higgs doublet <I >2 —)■ —$2 44,45 . Without loss of generality, 


one can then £x the sign of Ae; by convention, we choose Ag ^ 0 so that 0 ^ o; ^ 7r/2. 


2.2 Yukawa sector 

In the A2HDM, the interactions of the physical scalar helds with fermions are described by 
Cy = ~ H+\u[<;iVMiPR-^„lYllV PL]d + 

[fM,PRf] + h.c., (2.10) 


where Pl,r = (1 T 75)/2 are the chirality projectors, Mf=u,d,i represent the diagonal fermion 
mass matrices, V is the Cabibbo-Kobayashi-Maskawa (CKM) matrix, and 

Vd} — "^*1 + ("^*2 + i Tlis) ^d,i , 1/u* = Pil + {Tli2 — i Tlis) • (2-11) 


The parameters (/ = m, d, 1) are family-universal complex quantities which introduce new 
sources of CP violation beyond the CKM matrix. For particular real values of these parameters, 
indicated in table one recovers all different versions of the 2HDM with NFC. 

The Yukawa alignment condition is not stable against quantum corrections [7,13,48 . How¬ 
ever, the Favour symmetries of the A2HDM constrain tightly the possible FCNC effects, keeping 
them well below present experimental bounds [7 ,49-52 
duced at one-loop take the form 


The only FCNC local structures in- 


C 


C 


FCNC — 


47^2^3 


(1 + ^u'^d) ^ i -h iPjz){^d - ‘^u) 


- (71,2 - ^Pj3)(^d - 0 k VM^M^M^un 


dLV^M^MlVM^dR 
+ h.c., 


( 2 . 12 ) 
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Table 1: Two-Higgs-doublet models with natural flavour conservation. 


Model 

9d 


Q 

Type I 

cot fl 

cot (3 

cot (3 

Type H 

— tan /3 

cot (3 

— tan (3 

Type X (lepton specific) 

cot (3 

cot [3 

— tan (3 

Type Y (flipped) 

— tan (3 

cot (3 

cot f3 


which vanishes exactly for the 2HDMs with NFC. In general, the size of the indnced flavonr- 
changing interactions is controlled by three powers of qnark masses and the GIM mechanism. 
The renormalization of the conpling constant C is determined, nsing dimensional regnlar- 


ization, to be 17 


1 r 

C = CM + - I -p—^ +1E- 


(2.13) 


where yg ~ 0.577 is the Euler constant and /i is an arbitrary renormalization mass scale. The 
renormalized coupling satishes 


Ch(a) = Ch(Ao) - ln(A/Ao) ■ 


(2.14) 


Assuming Yukawa alignment to be exact at a given energy scale Aa, so that Cr{Aa_) = 0, 
implies that = ln(AA//i). 


3 Flavour-changing top decays 

The flavour-changing top decays t ^ ch and t ^ cV (V = 7 , Z) occur firstly at the one-loop 
level in the SM. The decay rate for these processes is not only suppressed by the loop factor, 
but receives in addition a strong CKM and GIM suppression Here we focus on 

hnal states with the charm quark because Br(t —)■ ■uX)/Br(f —)■ cX) ~ |I4b/IA6p ~ 7 x 10“^ 
in the SM, as well as in the A2HDM. Fixing the Higgs mass at — 125 GeV, one obtains 
the SM branching ratios: Br(t —)■ ch) ~ 0(10“^®), Br(t —)■ cy) ~ 0(10“^'^) and Br(t —)■ cZ) ~ 
0(10“^^). Within the A2HDM these decay rates can be enhanced due to additional charged 


Higgs contributions at the loop level. For t —>■ ap^ decays, the counter-term piece in eq. (2.12) 


would also contribute. 

In the SM, the dominant decay mode of the top quark is the unsuppressed two-body decay 
t M with F(t —)■ W^h)/mt ~ 1%. To compute the relevant branching ratios we take 


Br(t^cX) = 


T{t cX) 

rtot(^) 


(3.1) 


To a very good approximation, Ftot(^) — F(t —)■ W^b) holds also in the A2HDM, except when 
Mh± < rrit — rrih] in this case, the additional decay mode t — )■ H~^b must be taken into account. 
The partial decay widths for t —)■ W~^b and t M H^b are calculated at leading order: 


T{t^W^b) = 


9^\Vt 


tb 


6471 




Ml-) 


m. +mu + 


2 




(3.2) 
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r{t^H+b) = 


\Vt 


tb 


16 TT 




{ml + ml-Ml±) {ml\qd? + ml\q^\^) 


Amlml Re 


(3.3) 


Here A(a;, ?/, z) = + if + — 2{xy + xz + yz), and g is the SU(2)^ gauge coupling constant. 

Due to the smallness of rud^s and the unitarity of the CKM matrix, the t —)■ and t —)■ cV 

decay amplitudes turn out to be very sensitive to the bottom quark mass 33 . The most 


adequate choice for the internal quark masses is the running MS quark mass evaluated at a 
typical scale 0{mt) [^; we follow this prescription. The external quark masses are taken as 
the on-shell pole masses. 

In section phi] we describe the calculation of the t —)■ c(p° decay amplitudes. The decays t —)■ 


cV iy = 7 , Z) are discussed in section 3.2 Explicit analytical results for these processes within 


the A2HDM are collected in the appendices. All our results are presented in the Feynman gauge, 
corresponding to .^ = 1 in the gauge. We have however checked the gauge independence 
of our results by additionally performing all calculations in the unitary gauge. We have also 
checked analytically that our results reproduce the corresponding SM predictions [^|2^[^. 


3.1 t decays 

The total amplitude for t{pt) —t c{pc) decays can be written as Atot = A + A^t- Here A 

collects all the contributions arising from one-loop diagrams and can be parametrized as 27,29 

18 18 

.4 = 51 -4" = E E I’ll «(pt). (3.4) 


n=l 


n=l q=d,s,b 


Diagrams contributing to this process in the Feynman gauge are shown in figures [T] and The 
total loop amplitude is ultraviolet divergent, except in the particular cases of 2HDMs with 
NFC where a finite result is obtained. As expected, the divergence gets reabsorbed into the 


renormalization of the counter-term coupling C in eq. (2.12). The tree-level contribution of the 
counter-term Lagrangian takes the form 


4o. = 55 v;,r4 “(p.) 

q=d,s,b 


Oict Pr + /^ct Pl 


u{pt), 


(3.5) 


with 


ig^rrit 


«ct = - C 32^2]vf 3 + ^*u^d)y - y*{Pj 2 - iPjs) 


ig^rric 


y = - c ~ ^u)iPj 2 + iPjs) rUg 


(3.6) 


The partial decay width can be written as 


cv7°) = 


32 vrm? 


(|q!|^ -t- + ml - M\) 2mcmt {a*(3 + (3*a) 


(3.7) 
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Figure 1: Penguin diagrams contributing to t ^ in the Feynman gauge. 





Figure 2: Self-energy diagrams contributing to t ^ ap^ in the Feynman gauge. 


Here 


/ 18 

q=d,s,b \n=l 


18 


q=d,s,b 


ct 


(3.8) 


. n=l 


The contributions to the amplitude from each diagram, encoded in the coefficients and 
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are collected in appendix [b] 

3.2 t —cV (V = 7 , Z) decays 

The decays t{pt) —)■ c{pc) V{pv) can only proceed at the loop level (there are no tree-level 
counter-terms in this case). The decay amplitude can be parametrized as 


14 


14 


A = K,Kiu(p. 

n=l q=d,s,b 


n=l 




a) Py + ^2 Pi + as 7' 

p(^ + p'i + 7^1 Pnluipt) e*(pv), 


Pl 


(3.9) 


where is the polarization vector of the gauge boson V. The one-loop diagrams contributing 
to this process in the Feynman gauge are shown in figures and The total amplitude is of 
course ultraviolet hnite. 






( 3 ) 


( 4 ) 




( 6 ) 


( 7 ) 


( 8 ) 


Figure 3: Penguin diagrams contributing to t ^ cV {V = ^iZ) in the Feynman gauge. 

The partial decay widths for t ^ cV (17 = 7 , Z) decays are given, in the limit rric = 0, by 
r{t -^cj) = - [|a2|^ I&2P] - 2 [lasl^ Re(ata2) +m‘fRe{blb 2 ) 


2mtRe[a3(6i-h 62)] + 2 mtRe[bl{ai + a 2 )] >, 


(3.10) 


and 


F(t^cZ) = 


m; 


m: 


212 


128 vrm? 


3^2 1 ~ + l^2p] + 4 {nil 2M|) [|a3p I 63 P] 


-F 4mt {mf - M|) [Re(a 263 ) Re(a362)] 


(3.11) 
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Figure 4: Self-energy diagrams contributing to t ^ cV (V = 'y,Z) in the Feynman gauge. 


Here we have defined 


14 


14 


q=d,s,b \n=l 


(n) 


i>. = E E*.' 

q=d,s,b \n=l 


(n) 


(3.12) 


with i = 1,2,3. For the numerical analysis we always keep finite charm mass effects into 
account. The contributions to the amplitude from each diagram, encoded in the coefficients 
and b\^\ are collected in appendix C 


4 Discussion 

Relevant inputs for the evaluation of the flavour-changing top decay rates are listed in table 
We assume that the combined measurement of the top-quark mass by the Tevatron and the 
LHC corresponds to the pole mass, but we increase its systematic error by 1 GeV to account 
for the intrinsic ambiguity in the top-quark mass definition. The bottom- and charm-quark 
masses quoted in the table are MS running masses at the quark-mass scale, i.e., fnq{fnq), while 
the light-quark ones are MS running masses renormalized at a scale of 2 GeV. To fix the 
needed entries of the CKM matrix we use inputs with a minimal sensitivity to new physics 
contributions. Our SM predictions for the processes considered are presented in table 

The ATLAS and CMS collaborations have searched for flavour-changing decays of the top 
quark. The ATLAS collaboration sets the bound Br(t —)■ qZ) < 0.73% at the 95% confidence 
level (CL), with 2.1 fb“^ of data at ^/s = 7 TeV [^, where the q in the final state denotes 
a sum over q = u,c. The CMS collaboration has set a better limit, Br(t —)■ qZ) < 0.05%, 
with 24.7 fb“^ of data at ^/s = 7&8 TeV [^. The strongest current bound on f —)■ cy decay 
has been obtained by the CMS collaboration, Br(f —)■ cy) < 0.182%, using 19.1 fb“^ of data 
at ^/s = 8 TeV 39 The ATLAS collaboration sets the limit Br(f —)■ qh) < 0.79%, with 


^The CMS limit quoted here on Br(t —>■ cy) is actually derived from a search for the anomalous single top 
quark production in association with a photon in proton-proton collisions, following an effective Lagrangian 


approach with the assumption of vanishing contributions from both tqg and tuy interactions 


39 











Table 2: Relevant inputs for the evaluation of the decay rates. 


Input 

Value 

Gomment 

Input 

Value 

Gomment 

Gf 

Mw 

Mz 

7ckm 

IK.I 

iKfcl 

IKbl 

1.1663787(6) x 10“^ GeV”^ 
80.385(15) GeV 
91.1876(21) GeV 
(73.2+®;;))° 

0.2247(7) 

3.42(15) X 10-3 
42.21(78) X 10-3 

Ref. 

Ref. 

Ref. 

Ref. 

Ref. 

Ref. 

Refs, f 

54 

54 

54 

55 

56 

5*6 

7- 

59 

Mh 

rrit 

mb 

rric 

rris 

md 

Ois{Mz) 

125.14(24) GeV 
173.34(1.76) GeV 
4.18(3) GeV 
1.275(25) GeV 
93.8(2.4) X 10-3 GeV 
4.68(16) X 10-3 GeV 
0.1185(6) 

Refs. 

Refs. 

Ref. 

Ref. 

Ref. 

Ref. 

Ref. 

30, 

32 

54 

p 

56 

56 

54 

61 

64 


Table 3: SM predictions for flavour-changing top decays. 


Observable 

SM prediction 

Br(f —)■ cy) 
Br(f —)■ cZ) 
Br(f —)■ ch) 

(4.31 ±0.24) X 10-^^ 
(1.03 ±0.06) X 10-34 
(3.00 ±0.17) X 10-33 


25 fb ^ of data at y/s = 7&8 TeV 37 . A slightly stronger limit, Br(t —)■ qh) < 0.56%, has 
been obtained by the CMS collaboration, nsing 19.5 fb“^ of data at y/s = 8 TeV (^. Fntnre 


42,43 


prospects for these processes at the high Inminosity LHC have been discnssed in refs. 

One expects to improve the limits to the 10“^ level for Br(t —)• cV) {V = 'y,Z), while for 
Br(t —)■ ch) it would be possible to reach the 10“^ — 10“® level. 

For the phenomenological discussion we shall focus on the CP-conserving A2HDM, which 
contains 12 free real parameters: fi 2 , A*, {k = 1,...,7), the three alignment constants 
(/ = u,d,l) and the counter-term coupling CR^fi). Physical amplitudes are independent of 
the renormalization scale fx, due to eq. (2.14); in the following, we choose /r = Mw- Some 


of the parameters of the scalar potential can be traded by the physical scalar masses and the 
mixing angle a. The following relations 




T cos^ 


a -f sin^ 


a] , A4 = — \M? sin^ a -|- Mjr cos^ a -f M\ 

J yZ L 


■ 2M|±] 


As = \mI sin^ a -\- Mjr cos^ a — M^l , Xq = — Ml) cos a sin a , (4.1) 

2 v^ 


together with eq. (2.4), allow us to work with a set of parameters more closely related to 
physical quantities: 


Scalar sector: M^, Mr., M^, Mh±, cos a, A 2 , A 3 , Ay , 


Yukawa sector: Cr{Mw) ■ (4.2) 

Not ah the parameters are relevant for the processes we are concerned about. The decays 
t —)■ are only sensitive to {M^o, Mh±, cosd, A 3 , Ay, Cr{Mw)}. The transition 


9 
























































amplitude does not depend on the other neutral scalar masses M^o {i ^ j), as can be seen 
explicitly from eq. (B.2). There is also no dependence on the coupling A 2 ; the associated term 
in the scalar potential (<h 2 $ 2 )^ does not generate the needed cubic vertices because 


$2 has no vacuum expectation value (see eq. (2.1)). The decays t ^ cV {V = 7 , 2 ’), on the 


other hand, depend only on Mh± and the alignment parameters <;u,d- All the relevant cubic 
vertices are fixed in this case by the gauge symmetry and do not depend on free parameters of 
the scalar potential. 

We assume in the following that the 125 GeV Higgs boson corresponds to the lightest CP- 
even state h] he., we fix Mh — 125 GeV. The LHC data imply that it couples to the massive 
gauge vector bosons with a SM-like strength so that cosd ~ 1. We are interested in how 
large the enhancements of the flavour-changing top decay rates can be, compared with the 
SM predictions, focusing on the 125 GeV Higgs boson in the case of t —)■ c(p° transitions. To 
address this question, we analyze the parameter space of the A2HDM, subject to the following 
assumptions and constraints: 


The LHC and Tevatron Higgs data imply that cos a 
(/ = u,d,l) 1^,25. We work in the limit cosd = 


> 0.9 ( 68 % CL) and ||/J| ~ 1 
1 so that no constraints on the 


alignment parameters are obtained from the 125 GeV Higgs data 21,25 


We take into account constraints in the — (;d plane derived from the measurement of 


Br(i? —)■ Xs 7 ) 13,14 


We restrict the alignment parameter 


Z ^ bb decay and 


^ mixings 


< 2 , in order to satisfy the constraints from 
The parameters (;d,i are much less constrained 


13 


phenomenologically; we take < 50 as in ref. 14 


The four LEP collaborations, ALEPH, DELPHI, L3 and OPAL, have searched for pair- 
produced charged Higgs bosons in the framework of 2HDMs, excluding Mh± 80 GeV 


(95% CL) under the assumption that decays dominantly into fermions 65 


Searches for a light charged Higgs via the decay t —)■ H~^b performed by the ATLAS 66l67 
and CMS 


68,69 


collaborations, together with the limits on a charged Higgs from the 
Tevatron are taken into account. These direct searches give an upper bound on the 
Yukawa combination which, although being weaker than the one from Br(H —)■ 

Xs'f), basically exclude one of the two possible strips allowed by the latter (25). 


We consider the perturbativity bound on the quartic scalar couplings |A 3 ^ 7 | < dvr |^ . 
Additionally, the loop-induced decay h —)■ 77 is sensitive to A 3 and A 7 through the charged 
Higgs contribution to this process 21,25 . We take into account the latest measurements 
of the Higgs signal strengths in the h —)• 77 channel by ATLAS 71 and CMS 


In the limit cosd = 1, the decay rate for t ^ ch does not depend on Cr{Mw) and A 7 . 
Explicit expressions for all the relevant cubic Higgs couplings are provided in appendix In 
particular, for cosd = 1 we have A^+^_ = A 3 . The measured Higgs signal strengths by ATLAS 
and CMS in the di-photon channel are then only sensitive to A 3 and Mh±. Since in this case 
the Higgs production cross-section is the same as in the SM, one can write the Higgs signal 
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strength in the di-photon channel as I^I^: 


fi. 


'77 


a{pp ^h)x Br(h ^ 27 ) ^ V 

(t(pp-)■ h)sM X Br(h 27 )sm ^ 


where encodes the charged Higgs contribution to h —)■ 27 and is given by 


(4.3) 


/^h _ 

o^± - 




Here 


X 


A{x) = -X- ^/(a:), 


^H+H- x4.(a;/f±). 


f{x) = —4arcsin^(l/\/a;), 


(4.4) 


(4.5) 


with xh± = We require that the Higgs signal strength in eq. (4.3 

range of the experimental measurements. The latest results by ATLAS |71 


lies within the 2 cr 


/i: 


77 


_ 1 1 7+0.28 

~ ' - 0 . 26 > 

and by CMS |^: = 1.12 ± 0.24, are consistent with the SM. 

Performing a scan over Q, A 3 }, subject to the restrictions specihed above, while 

hxing the charged Higgs mass to benchmark values, we obtain the upper bounds on Br(f —)■ cV) 
{V = 7 , Z) and Br(t —)■ ch) shown in table In the window 90 GeV < Mh± < 150 GeV the 
alignment parameter qj. is constrained to be small by the direct charged Higgs searches at the 
LHC via top decays, |<^rf| < 10, implying a very strong suppression on the decay rates. For 
Mh± < 90 GeV a weaker bound on |<^d| is obtained by a combination of LHC and Tevatron 
limits, |<^d| < 25. For Mh± > 150 GeV the largest decay rates for these processes are obtained 
for |<^„| < 1 and ~ 50. The upper bounds obtained for Br(f —>■ cV) put these processes well 
beyond the reach of the high luminosity LHC, within the A2HDM (4^. Similar conclusions 
were obtained in refs. 


34,35 within the framework of 2HDMs with NFC. 


The decay rate for t ^ ch can receive on the other hand much larger enhancements, due to 
the intermediate charged Higgs contribution involving the cubic Higgs coupling The 

maximum values for Br(f —)■ ch) are obtained when the cubic scalar coupling A^+j^_ saturates 
either the h —>■ 27 limits or the perturbativity bound. Diagram 3 in hgure dominates the 
corresponding decay amplitude in this case. The contribution from this diagram to the decay 
amplitude is proportional to qu^d^^+H- table While the product quqd is 

constrained to be small in magnitude by Bt{B —)• Xg'y), the term proportional to becomes 
greatly enhanced for large |<^d| values. Such large values of |<^d| can be obtained outside the 
window 90 GeV < Mh± < 160 GeV since the limits from direct charged Higgs searches via top 
decays at the LHG are avoided. 

Analyses of f —)■ ch decay within the type H 2HDM, prior to the Higgs discovery, have 
found that a light charged Higgs can enhance considerably the associated decay rate in this 
model, for large values of tan/3 and the cubic Higgs coupling and even reach the level 

of expected sensitivity at the high luminosity LHG: Br(f —)■ ch) ~ 10“® 34,35 . Such behavior 
is compatible with our hndings, given that in the limit qd = ~ — tan/3 we recover the 

Yukawa couplings of the type H 2HDM. However, we hnd that current measurements of the 
125 GeV Higgs properties play an important role when evaluating possible enhancements of 
Br(f —)■ ch). In particular, measurements of the Higgs signal strengths in the di-photon channel 
restrict the allowed size of the cubic Higgs coupling X^+jj- for a light charged Higgs. This in 
turn implies that the allowed enhancements of Br(f —)• ch) cannot be as large as previously 
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Table 4: Upper bounds for Br(i —>■ cV) (V = ^,Z) and Br(t —>■ ch) in the CP-conserving A2HDM. 


Mh± [GeV] 

Br(f —)■ cy) 

Br(f —)■ cZ) 

Br(f —)■ ch) 

100 

< 2 X 10-^2 

< 2 X 10-^3 

< 6 X 10-1^ 

200 

< 10-^° 

< 3 X 10-11 

< 3 X 10-« 

300 

< 10-11 

< 5 X 10-12 

< 2 X 10-® 

400 

< 2 X 10-^2 

< 2 X 10-12 

< 5 X 10-1* 

500 

< 10-12 

< 10-12 

< 2 X 10-9 

Exp. limit 

< 1.8 X 10-^ 39 

< 5 X 10-4 38 

< 5.6 X 10-*i 40 


speculated. Taking into account the measurements of the 125 GeV Higgs properties, searches 
for a light charged Higgs via top decays, and the flavour constraints specihed earlier, we hnd 
that the decay rate for t ^ ch lies beyond the reach of the high luminosity LHC in 2HDMs 
without tree-level FCNCs. Under the constraints considered the largest decay rate is obtained 
for Mh± being slightly below 90 GeV, Br(f —>■ ch) < 2 x 10“^. 

It is necessary to discuss the robustness of the previous statement. If small deviations 
from the limit cosd = 1 are considered, the LHC Higgs data gives rise to strong bounds on 
the magnitude of the alignment parameters. Since \yj\ = |cos(5 + <^/sind| (/ = u,d,l) is 


constrained to be close to one, one obtains |<^/| < 0{1) when cosd < 1 This implies 

in particular that |<^d| should be small and large enhancements of Br(f —)■ ch) are not possible. 
Allowing for CP violation would not led to any signihcant enhancement either, given the strong 
constraints on CP-violating couplings derived from electric dipole moment experiments [l6| . 

We turn now to discuss the role of the direct counter-term contribution to t ^ ch decay, 
which is not present in 2HDMs with NFC. In the limit cosd = 1 this contribution vanishes 
because of the orthogonality of IZ. The LHC data imply that cos a is very close to one so that 
the counter-term contribution to the flavour-changing t ^ ch decay will be suppressed by a 
small factor sind at the amplitude level. Furthermore, the characteristic flavour structure of 


the A2HDM counter-term (2.12) implies a strong suppression of its effects, due to the explicit 


powers of quark masses and the unitarity of the quark mixing matrix [^. Neglecting the loop 
contribution (at y = M]y), 


Br(t y ch)^j-ee 


(1 - Mi/m 




2sm^9wM^ (1 - (1 + 2M^/m2) 

2 X 10“^^ sin^dlUrfp, 


sin^ a\Ed\ 


(4.6) 


where 


Ed — Cfi{Mw) (1 + ‘iu‘id) {‘id — ‘iu) ■ 


(4.7) 


The size of Ed is constrained experimentally by the measured amount of mixing between the neu¬ 


tral Bi meson and its antiparticle, which receives also contributions from the Lagrangian (2.12), 


mediated by the three neutral scalars ip/ = {h,H,A}. One hnds that this process allows for 
\Ed\ ~ even when the masses of the neutral scalars are of 0(100 GeV) 1^, but this is 
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far too small to generate any observable signal in t —)■ c/ij^ 

It is important to analyze also the impact of the perturbativity bound on the results ob¬ 
tained. A different upper limit on the relevant cubic coupling \ has been considered in 

ref. 21 . The charged Higgs gives the following hnite correction to the hH^H~ vertex, at the 


one-loop level; 


i^H+H-)eS — ^H+H- 


1 + 


167r2M^± 




Ml 

Ml. 


= ^H+H- (1 + 


(4.8) 


where ^ ^ 

Z{X) = f dy f dz [{y + zf + X{l-y-z-yz)]~^. (4.9) 

Jo Jo 

Large values of |A^+^_ | could make this loop correction comparable to the leading-order result, 
which would cast doubts on the perturbative expansion. We therefore allow this correction to 
be at most 50% (A ^ 0.5). In figure]^ we show the region allowed at 2cr by the measurement of 
the Higgs signal strengths in the di-photon channel, together with the bounds extracted with 
the perturbativity limits IA 3 I ^ dvr and A ^ 0.5. 



Figure 5: Allowed region by measurements of the Higgs signal strengths in the di-photon ehannel 
(blue-meshed) together with the perturbativity limits IA 3 I ^ dvr (light gray) or A ^ 0.5 (dark gray). See 
text for details. 


The constraints from h —)■ 77 give rise to a large allowed region, centered around A 3 = 0, 
whose width increases for higher values of Mh±. In this area the h —)■ 77 decay amplitude is 
dominated by the W-boson and top-quark loop contributions, as in the SM; the charged Higgs 
contribution remains subdominant. For light charged Higgs masses a small disjoint allowed 
region appears with A 3 > 6 . In this small area the charged Higgs contribution dominates over 


the PF-boson and top-quark loops and flips the sign of the amplitude, (1 — 0.15 Cj^i) ~ —1, 
giving a SM-like Higgs signal strength (see eq. (4.3)). In principle there is no reason to expect 
such an accidental tuning of the charged Higgs contribution to occur. This separate small 


2 Assuming Yukawa alignment to hold at the high-energy scale A^, we have Cr{Mw) = \n{AA/Mw). 
Therefore, Cij(A'/w)/(47r2) < 1 for A^ < 10^® GeV. 
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region is therefore not to be seen as very realistic. It is possible to argue that such region 
brings problems to the perturbative expansion. The perturbativity limit A ^ 0.5, being more 
stringent for light charged Higgs masses, excludes this small region. For a light charged Higgs 
the maximum values of Br(f —)■ ch) are obtained precisely in this separate region, where the 
value for | reaches its maximum allowed value. The limits on Br(t —)■ ch) would therefore 

be even stronger for a light charged Higgs if the perturbativity limit A ^ 0.5 is taken into 
account. Once this perturbativity limit is considered, we get the limit Br(t —)■ ch) ^ 6 x 10“®. 

Recent works have studied the possibility to look for flavour-changing top-quark anomalous 
interactions via production processes. In ref. 72 a fully gauge-invariant effective-held-theory 


approach was adopted for parametrizing the top-quark FCNC interactions, while in ref. 73 


it was assumed that the Higgs boson posses tree-level flavour-changing couplings with the top 
quark. It was pointed out in these works that these top-quark flavour-changing effects can 
often be probed with better sensitivities in production processes than via the top-quark decays. 
Whether this is also the case within the A2HDM with a light charged scalar deserves a detailed 
analysis but lies beyond the scope of the present work. 

5 Conclusions 


We have performed a complete one-loop calculation of flavour-changing top decays {t —)■ cy, 
t —)■ cZ, t —)■ c</ 9 °), within the A2HDM. Here = {h,H,A} represents any of the neutral 
scalar mass eigenstates. Our results agree with the available SM results in the literature when 
the corresponding limit is taken 27,29 . We have also checked the gauge independence of our 


results by carrying out the calculation in the Feynman and unitary gauges. Explicit analytical 
expressions in the Feynman gauge are provided in appendices|^and[^ The results are presented 
in the limit rric = 0 , in order to avoid lengthy expressions; however, in our numerical analyses 
we have always used the exact expressions. 

The SM predictions for these transitions are given in table They are orders of magnitude 
too small to be accessible even at the high-luminosity phase of the LHC. We have investigated 
whether signihcant enhancements of the branching ratios could be possible within the A2HDM. 
Assuming that the 125 GeV Higgs-like boson corresponds to the lightest CP-even state h of the 
CP-conserving A2HDM, we have discussed the impact of the relevant model parameters on the 
decay rates. We have taken into account the constraints from flavour experiments as well as 
the measurements of the Higgs-boson properties at the LHC. Upper bounds obtained for these 
rare top decays within the A2HDM are listed in table for benchmark values of the charged 
Higgs mass. 

While sizeable enhancements are indeed possible, compared with the SM predictions, we 
hnd that the decay rates for t —>■ cV {V = 7 , Z) remain well below the expected sensitivity 
levels at the high luminosity LHC, across all of the parameter space considered. As long as the 
charged Higgs is relatively light, the decay rate for t ^ ch receives much larger enhancements 
for large values of and the cubic Higgs coupling The LHC measurements of the Higgs 

signal strengths in the di-photon channel are found to play a very important role in estimating 
the maximum allowed values for Br(f —)■ ch). The charged Higgs also contributes at the loop 
level to the decay h —)■ 77 and, for a light charged Higgs, a large cubic coupling would 

led to large deviations of the Higgs signal strengths in the di-photon channel. We hnd that 
Br(t —)• ch) lies also beyond the reach of the high luminosity LHC, once the constraints from 
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the Higgs data are taken into acconnt. 
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A Loop functions 

Dimensional regularization is used in our calculations. The scalar loop functions appearing are 


given by 74 


Ao(Mi) 

Co(/,S,Mi,M2,M3) 


= y d^k 

= j d^k 

= f d^k 


- Ml ’ 


(P -M 2 )p + /)2 -M|] ’ 

1 


Here 


(F - Ml)[{k + If - Ml][{k + / + s)2 - M|] 
d^k 


d‘’h = 


(2w: 


D 


(A.l) 


(A.2) 


is the integration measure with e = 4 — D, and is the SU(2)^ gauge coupling constant 

in D dimensions. Vector and tensor integrals are reduced to the scalar loop integrals via the 
Passarino-Veltman method 75 . Following the notation of refs. [^[^, we have 

Bi = ^ [Ao(Mi) — Aq{M2) — SiRo] ) 

Cq = Bq{s, M2, M3) + M^Cq , 

where Si = P + M^ — M|. The other relevant loop functions are given by 


(A.3) 


Cii 

C 12 

C 21 

C 23 


= Y 


= Y 


Bq{1 -t- s, Ml, M3) — Bq{s, M2, M3) — SiCo 
Bq( 1, Ml, M2) — Bq {1 -f s. Ml, M3) — S2CQ 

Bi{l -t- s. Ml, M3) + Bii{s, M2, M3) — siCii — 2(724 
Bi{l, Ml, M2) -Bi{l + s. Ml, M3) - S2C11 


(A.4) 
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and 


C 22 — 


-I ■ s 


Bi{l + s, Ml, M 3 ) — -Bi(s, M 2 , M 3 ) — S 1 C 12 


C 2 A — 


+ 1“^ ~Bi[l + S, Ml, M3) — 52^*12 ~ 2C24 

Bo{s, M2, M3) + 2 Mi Co + siCii + S2C12 


2(C)-2) 

Here we have defined S 2 = + 2/ ■ s + M| — M| and 

1 


(A.5) 


Y = 


2 [Ps^-ii-sy] 


-l-s 


-l-s f 


(A. 6 ) 


B Decay amplitude for t 


C<p° 


We parametrize the one-loop contribution to the t ^ decay amplitude as indicated in 


eq. (3.4). In tablej^we give the analytical expressions for the coefficients and obtained 
from the 18 Feynman diagrams in figures [T] and For simplicity, we only give the results in 
the limit me = 0 , although we have used the exact expressions, including finite charm masses, 
in our numerical results. In this limit, all coefficients = 0 (n = 1,..., 18), while = 0 
for n = 15,16,17,18. 

We have defined the combination 


Id' = B-Y + ^12 Re(<^d) - 77j3 Im(<^d). 


(B.l) 


The matrix 77 determines the neutral Higgs boson states in terms of the neutral components 


of the scalar doublets in the Higgs basis (see sec. The parameters appearing in table 
have been defined in eq. (2.11). The relevant cubic couplings in this case read 17 




\A — — 7? 

w+w- G+w- ’ 

^H+W- ^ ”^12 — i B-js , 


A 


H+H- ~ -^3 B-ji + Af TZj2 — X\TZj^ , 


M\ 


^G+G- ~ ^"^3 ^72 ~ Af; 77,3 — 


L-jl 




6 '^3‘i 


'6 'M3 


77ji, 


0 1 i M^o-M^i 

^H+G- ~ -^6 + 2 (-^4 + 2A5) 77j2 — - (A4 — 2A5) 77j3 = - 


j2 ^ ”^ 73 ) 


(B.2) 


R I 

Here A^ ’ denote the real and imaginary parts of Xk respectively. 
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C Decay amplitude for t ^ cV 

Following the notation of eq. ( |3.9[ ), all non-vanishing contributions to the t —?• cV (V = 7 , Z) 
decay amplitude have been given, for rUc = 0 , in tables id and|^ Here we have dehned 

gvq = Avq + Bvq''f5, (C-l) 

and 

Sq = Avq + Bvq , Pq = Ayq — Byq , (C. 2 ) 

with q = (m, d). Values for the relevant constants are given in tables]^ andThe weak mixing 
angle is determined by e = g sin 9w, Mw = gv/2 and My/ = Mz cos 9w We use the notations: 
sw = sin^^y and cw = cos9w 
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Table 5: Amplitude for t ^ ccpj in the limit me = 0. 
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(n) 
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Table 6: Amplitude for t ^ cV {V = Z) in the limit me = 0; Coefficients 


n 

(n) 

03 

argument 

1 

• 2 2 

{PV, -pt,mq,mq,MH±) 


— — 2?*i?yrfmi] (Cii — C 22 ) + 2<^^5rfC'24| 


2 

• 9 2 

ig r ~ . 0 Ox 

2M2 I ‘S'dC'o + Pd(mg mj2)C'o 

{Pv, -pt,mq,mq,Mw) 


+ [2-Bydm^ — MySd\ (Cn — ( 712 ) + 2SdC2^ 



,.2 9vh+h-^24 


3 

{Pv, -Pt, Mh± , Mh± , mq) 

4 

g‘^mq 

I\j2 9VG+G-C24 

{pv, -Pt, Mw, Mw, mq) 

5 

. 2 

- 2)PdCo -{D- 2)mpdCo + ^4(11 - 2)M^Cn 

{Pv, -pt,mq,mq,Mw) 


- {{D - 4)M2 + 2ml) C 12 - 2(11 - 2 )^ 24 ] } 


6 

2^^ S'V'G+W-jo^gC'o + mac'll <^ 12 )! 

{pv, -Pu Mw, Mw, mq) 

7 

g'^mf 

{pv, -Pt, Mw, Mw, mq) 

8 

^9vw+w- 1 “2(7o + 2{ml — My)Co 
+{ml - 2M2 )Cii + mlCu - 2{D - 2)C724} 

{pv, -Pt, Mw, Mw, mq) 

9 

,g‘^ml 

{pv - Pt,MH±,mq) 

10 

g'^mf 

i ^ PuBq 

2M4 

{pv -Pt,Mw,mq) 

12 


{-Pt,MH±,mq) 

13 

g'^ml 

i PuBi 

2M4 

{-Pt,Mw,mq) 

14 

'-{D-2)g^P^{Bo + Bi) 

{-Pt,Mw,mq) 
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Table 7: Amplitude for t ^ cV (F = 7, Z) in the limit me = 0; Coefficients 


n 


argument 

1 

Wo + - CuPd)Cn - CuSdiCu - C 21 + C 23 )} 

{Pv,-Pt,mq,mq,MH±) 

2 

-i {PdCo + SdiCu - C 21 + C 23 ) - 2BvdCii\ 

{pv,-Pt,mq,mq,Mw) 

3 

2m2 + 2?„)C'ii ?rf(C'i2 2(721 + 2(723)1 

{Pv , -Pt , Mh± , Mh± , mg) 

4 

-^-J^9vg+g- {Co + SCii - C 12 + 2 C 21 - 2 C 23 ) 

{Pv, -Pt, Mw, Mw,mq) 

5 

ig^mtPd[{D - 2)(7ii - {D - 4)(7i2 + {D-2) [C 21 - C23\ } 

{pv,-Pt,mq,mq,Mw) 

6 

M^3vg+w-{Co + Cu) 

{pv, -Pt, Mw, Mw,mq) 

8 

2 

^-^9vw+w-{^Co -DCii + {D + 2)Ci2 - 2{D - 2) [C 21 - C 23 ] } 

(pv, -Pt, Mw, Mw,mq) 


Table 8: Amplitude for t ^ cV (V = j, Z) in the limit me = 0; Coefficients 


n 

"2 

argument 

1 

9d{<,lPdCo + 5d«Cl2 + ?^C 22 - ?^C23) - 2BvdCuCll} 

{Pv,-Pt,mq,mq,MH±) 

2 

iq^mfmt r 'i 

m 2 {PdCo + SdiCi 2 + C 22 - C 23 ) - 2BvdCn} 

{pv,-Pt,mq,mq,Mw) 

3 

m 2 9 9VH+H-^d\^‘iuCl2 + 9a{C23 (722)| 

{pv, -Pt, Mh± , Mh± ,mq) 

4 

q^mimt / \ 

m 2 9VG+G- [C 12 C 22 + C23j 

{pv, -Pt, Mw, Mw,mq) 

5 

- ig^mtPd{2Ci2 -{D-2) [C 22 - C 23 ] } 

{pv,-Pt,mq,mq,Mw) 

6 

M^ 9vg+w- l<-"0 + +iij 

{pv, -Pt, Mw, Mw,mq) 

8 

9‘^'^t9vw+w- ~ Co — 2 Ci2 — (D — 2) [C 22 — < 7 * 23 ] 1 

{pv, -Pt, Mw, Mw,mq) 
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Table 9: Quark couplings with neutral vector bosons. 



Table 10: Cubic couplings of the neutral vector bosons. 


V 

gVH+H- 

9vg+g- 

9vg+w- 

9vw+w- 

7 

e 

e 

—eMw 

e 

Z 

e cot{29w) 

e cot{29w) 

gswMz 

gcw 
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